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Abstract 


Various candidate aircraft and spacecraft materials were analyzed and com- 
pared in a low-energy neutron environment using the Monte Carlo N-Particle 
(MCNP) transport code with an energy range up to 20 MeV. Some candidate 
materials have been tested in particle beams, and others seemed reasonable to 
analyze in this manner before deciding to test them. The two metal alloys ana- 
lyzed are actual materials being designed into or used in aircraft and spacecraft 
today. This analysis shows that hydrogen-bearing materials have the best shield- 
ing characteristics over the metal alloys. It also shows that neutrons above 1 MeV 
are reflected out of the face of the slab better by larger quantities of carbon in the 
material. If a low-energy absorber is added to the material, fewer neutrons are 
transmitted through the material. Future analyses should focus on combinations 
of scatter ers and absorbers to optimize these reaction channels and on the higher 
energy neutron component ( above 50 MeV). 


1. Introduction 

Ionizing radiation is a concern to the occupants 
and equipment of aircraft and spacecraft. The interior 
environment is modified from the ambient by the 
structure of the craft and its contents, which include 
the flight crew and passengers, if applicable. There- 
fore, the skin material and structure of the craft is an 
important component for determining the occupant 
and equipment radiation exposure along with its pri- 
mary structural role. 

Neutrons are generated when high-energy parti- 
cles interact with matter. Because radiation sources for 
aircraft include the particles generated in the atmo- 
sphere above the aircraft, more neutrons exist in this 
situation than in a spacecraft environment. This fact 
does not mean that neutrons in the spacecraft environ- 
ment are not important, just that a different modeling 
approach is needed because all the neutrons are gener- 
ated in the spacecraft skin. Therefore, the focus of this 
paper is toward the aircraft environment. 

The ambient radiation fields at most altitudes in 
the atmosphere are fairly weak, but the cumulative 
effect over the number of flight hours is where the 
concern lies. Also, as the flight altitudes increase for 
subsonic aircraft (47 000 to 55 000 ft) and especially 
for supersonic aircraft (60000 to 70000 ft), the ambi- 
ent radiation levels increase and the flight crew, pas- 
sengers, and electronic equipment become more 
susceptible to solar particle events. 


Flight crew members are currently limited to 
1000 flight hours per year, but passengers have no 
limits. To put this in perspective, current regulations 
for radiation workers allow nonpregnant personnel to 
receive 5 rem on an annual basis and the general pub- 
lic is limited to 0.1 rem in all situations. Recommenda- 
tions are being considered (refs. 1 and 2) to lower the 
limit to 2 rem for nonpregnant radiation workers. It 
must be noted, however, that U.S. flight crews are not 
currently considered radiation workers. Therefore, 
they have no direct radiation exposure limits or moni- 
toring requirements, just a limit of 1000 flight hours 
per year. 

The high altitude, supersonic aircraft environment 
is also similar to the radiation environment measured 
at the surface of Mars. The density altitude of these 
aircraft is 20 to 25 g/cnr, and the density altitude for 
the Martian surface is 20 to 25 g/cnr. The constituent 
difference in the two atmospheres (air for Earth and 
CCF for Mars) is of little relevance. Therefore, any 
information gathered about high-altitude aircraft mate- 
rials is applicable to materials needed for the habita- 
tion of Mars. 

A large portion of the high-altitude radiation envi- 
ronment consists of neutrons (ref. 3). These neutrons 
appear mainly in two energy regions: a lower energy 
region between 0.1 and 10 MeV and a higher energy 
region between 50 and 1000 MeV. The low-energy 
neutrons are from an evaporation process where a neu- 
tron is boiled off from the nucleus because it contains 
an excess of energy from its last interaction. The high- 
energy neutrons occur from high-energy collisions 



(ten's of GeV/nucleon) where the result of the interac- 
tion is a spray of nuclear fragments and particles. 

The neutron evaporation spectrum has the same 
characteristics as a fission spectrum minus the gamma 
rays. Therefore, a standard nuclear reactor engineering 
model with experimentally derived interaction proba- 
bilities can be used to computationally test various 
materials in a fission spectrum to detemiine their neu- 
tron radiation modification characteristics. These 
results can then be used initially to evaluate, for the 
lower energy neutrons, various material combinations 
which can be compared with experimental data when 
available. 

2. Material Selection and Simulation 

The Monte Carlo N-Particle (MCNP) (ref. 4) pro- 
gram is used to simulate the particle interactions based 
on a Monte Carlo solution technique to the Boltzmann 
transport equation. The geometric model is a finite, 
plane-parallel slab, 3 g/cm 2 thick at 300 K. The source 
is an isotropic planar source at the left face of the slab 
and is modeled as an unmoderated Cf neutron 
source. Nine sample materials, enumerated in table 1, 
are each modeled as a single homogeneous slab. These 
nine materials represent stmcture or shielding materi- 
als currently being tested or materials that might be 
useful for stmcture or shielding but have never been 
analyzed or tested. The MCNP ID column lists the 
cross section sets used, which are the latest evaluated 
nuclear data file B (revision) (ENDF/B) evaluations 
available. This analysis method allows a comparison 
of the different materials without actually trying to 
model a specific component for a specific task. Impor- 
tant reaction channels can be identified within a mate- 
rial and an optimal shield can be designed to either 
encourage or discourage a particular channel through 
material selection and placement. 

The aluminum alloy (No. 2024) is a standard light 
alloy being used in various aircraft and spacecraft 
structures. The titanium alloy (Ti6A14V) is being con- 
sidered for future high-altitude aircraft. The rest of the 
materials are a polymer-based carbon-fiber-reinforced 
epoxy composite (AS4/3502 graphite/epoxy compos- 
ite) and polyethylene materials with various additives. 
The boronated polyethylene is a standard neutron 
absorber used in the nuclear industry. The beryllium 


and lithium additives are being investigated for their 
scattering and absorption properties. 

The polymer-based composite material is being 
investigated for various reasons. First, these compos- 
ites can be used as stmcUiral materials like alum in um 
and titanium. Second, they are lighter and therefore 
require less fuel to get to high altitudes and orbit. 
Finally, the reaction channels for nuclear fragmenta- 
tion for the lighter elements in the shield are less 
numerous than for the heavier metal alloy elements. 
Therefore, the lighter materials should not create as 
many new particle fragments. 

3. Results and Discussion 

MCNP tallies were created to obtain the energy 
spectra for entrance and exit surface fluxes, total reac- 
tion rate, and elastic scattering and total absorption 
reaction rates in the slab material. Also, energy- 
integrated values of these tallies were created. With 
this information, a description of the differences 
among materials and why these differences occur can 
be formulated. 

The source distribution is a typical Cf 252 spectrum 
inherent in MCNP and is shown in figure 1. This 
source peaks at an energy of 2.5 MeV with its maxi- 
mum energy on the order of 25 MeV. The evaporation 
spectrum of neutrons at high altitudes mimics this 
spectrum. Therefore, this source distribution can be 
used to approximate the evaporation spectrum that the 
materials in table 1 would encounter at high altitudes. 

The energy spectra for the transmitted flux values 
are shown in figures 2 and 3 and are normalized by the 
Cf source flux from figure 1. From these spectra, 
the hydrogen-containing materials are seen to transmit 
fewer neutrons with energies above 100 keV than alu- 
minum and titanium alloys; however, they transmit 
many more neutrons below 100 keV than even the 
source does. Fortunately, the biological damage 
(energy available to transmit to tissues) associated 
with these lower energy neutrons is small in compari- 
son with those neutrons with energies above 100 keV. 

The hydrogen-bearing materials reflect more low- 
energy neutrons than the metal alloys as shown in fig- 
ures 4 and 5. Because the hydrogen nucleus (a proton) 
is approximately the same mass as the neutron, it is the 
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best neutron-scattering medium available and is the 
driving reaction channel for neutron reflection. An 
increase in the number of particles reflected out of a 
shield means the fewer particles that are available to 
affect tissues and cause biological damage. However, 
reflected neutrons from a vehicle in a space environ- 
ment now act as a neutron source for extravehicular 
activities. This fact must be considered when some- 
thing like the Space Station is being constructed by 
humans outside the vehicle. 

Of the polyethylene-based materials, the prue 
material has the lowest transmission flux above 
100 keV but the highest flux below 100 keV in both 
figures 2 and 3. This could be because the other mate- 
rials have neutron scattering and absorption materials 
incorporated into their matrices displacing the hydro- 
gen. The investigation of the reaction channels will 
clarify whether absorption or scattering (leakage) are 
causing these differences. 

From the energy-integrated flux values listed in 
table 2, it is apparent that the hydrogen-rich materials 
also transmit the fewest total number of neutrons 
through the shield even though this is not apparent 
from figures 2 and 3. Additives to the basic polyethyl- 
ene base have some effect on the integral shielding 
properties. For example, if an efficient neutron 
absorber is added, like B 10 , then the total number of 
transmitted neutrons is reduced by another 10 percent. 

The reflected flux is a measure of the total scatter- 
ing properties of the material. Again, the more light 
elements in the material, the larger the integrated 
reflected flux will be. It is worth noting that of the 
polyethylene-based materials, the B 1 0 material pro- 
duced the smallest integrated reflected flux, probably 
because the number of neutrons available for the scat- 
tering reaction channel is reduced because of the high 
absorption. The metal alloys produce the smallest inte- 
grated reflected flux with the composite material 
between the two groups. The metal alloys are large 
atoms and are inefficient scattering materials. The 
composite material contains more carbon than hydro- 
gen. Carbon is a good scatterer but not as efficient as 
hydrogen. 

The spectral infomration in figures 4 and 5 shows 
the expected results and an unexpected one. Again, the 
polyethylene-based materials scatter the most neutrons 
to lower energies and allow them to backscatter out of 


the material at those energies. However, for neutrons 
above 1 MeV, the aluminum alloy backscatters more 
neutrons than every other material except the compos- 
ite. The larger titanium atoms would be expected to 
backscatter more than aluminum or the composite at 
those energies. However, the neutron absorption of the 
titanium alloy is much larger above 1 MeV than that of 
the composite or aluminum (ref. 5). 

Figures 6 and 7 show the expected results for the 
total scattering reaction rate energy spectrum. At ener- 
gies above 1 MeV, the titanium alloy shows a peak 
above all other materials because of the relatively 
large density of the titanium alloy compared with the 
other materials. However, the scattering reaction rate 
for both metal alloys drops very steeply as the energy 
of the neutron decreases because of the size of these 
atoms and the inability of these atoms to reduce the 
energy of the neutron per interaction. For large atoms, 
the energy reduction per interaction is small; there- 
fore, the metal alloys do not scatter neutrons to the 
lower energies. 

The hydrogen-containing materials, however, 
spread the neutron population over all energies in the 
spectrum because of the large average energy loss per 
interaction. The additives introduced into the hydroge- 
nous material (boron, beryllium, and lithium) ulti- 
mately reduce the number of neutrons at the lower 
energies because of the reduced density of hydrogen in 
the materials and the absorption properties of the addi- 
tives. The material that spreads more neutrons over all 
energies is pure polyethylene. This material is also the 
one with the largest hydrogen density. The composite 
material has the lowest hydrogen density, and even 
though neutrons are spread over all energies, the num- 
ber of neutrons at the lower energies is an order of 
magnitude lower than for polyethylene. 

Figures 8 and 9 show the total absorption reaction 
rates for the entire slab, which is important because it 
can keep neutrons from penetrating the shield and con- 
tributing to the dose of the occupants. These figures 
show typical absorption reaction rates for these mate- 
rials. The metal alloys show many more resonance 
peaks over the spectrum than the hydrogenous materi- 
als. Also, the metal alloys do not scatter as many 
neutrons to lower energies; therefore the reaction rates 
decrease at the lower energies. Of the hydrogen- 
bearing materials, the materials with boron- 10 and 
lithium- 6 absorb many more neutrons than any of the 
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others. Also, these two materials show low-energy 
absorption peaks rather than a steady rise just because 
the flux increases at these energies. Therefore, these 
materials are demonstrating true absorption, instead of 
a flux-related absorption. 

The reaction rates for (n,2n) and (n,3n) were so 
small, that they have not been included in this assess- 
ment. These reaction channels become much more 
important when the energy of the neutron becomes 
much larger (greater than 50 MeV). 

Based on the foregoing analysis, pure polyethyl- 
ene has the largest flux below 100 keV because it is 
the largest scatterer and smallest absorber of all the 
materials considered. Therefore, pure polyethylene is 
not considered to be the best choice for a shielding 
material. Borated polyethylene appears preferable, but 
a combination of various materials might outperform 
it. Clearly, considerable work remains yet to be done. 

An interesting point to note is that the broad 
absorption resonance from boron (fig. 8) starting at 
about 100 keV does not overlap a pronounced peak 
from lithium at 300 keV (fig. 9). Also, the composite 
material is the best reflector of higher energy neutrons. 
These three materials could be combined together in a 
shield design to produce a better shield than any of 
them alone. Of course, the high-energy (greater than 
50 MeV) response of this new shield would have to be 
analyzed to make sure these material additives do not 
increase the dose from their respective components. 

4. Concluding Remarks 

The hydrogen-bearing materials reduce the num- 
ber of higher energy neutrons in the 500 keV to 


20 MeV energy range that penetrate the shield materi- 
als considered. In fact, the more hydrogen, the better 
the overall shielding characteristics. These results 
also show that if neutron absorbers are added to the 
hydrogen-bearing material, then fewer neutrons pene- 
trate the shield. If a polymeric composite material is 
added, then a better higher energy reflecting character- 
istic could be achieved. These results need to be 
extended and refined to design a shield material that 
can optimally shield for low-energy neutrons. This 
new shield then needs to be analyzed for the high- 
energy neutron component (greater than 50 MeV). 
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Table 1. Elemental Compositions of Materials Analyzed 


Material 

Density, 

Thickness, cm 

Constituents 

MCNP ID 

Atomic 


g/cm 3 

at 3 g/cm 2 

fraction 




Al-27 

13027.60c 

0.944469 




Li-7 

03007.60c 

0.035748 

2024 A1 alloy 

2.7 

1.321586 

Cu-63 

29063.60c 

0.011679 



Cu-65 

29065.60c 

0.005206 




Li-6 

03006.60c 

0.002899 




Ti-nat 

22000.60c 

0.876018 

Ti6A14V alloy 

4.46 

0.672646 

Al-27 

13027.60c 

0.069090 



V-nat 

23000.60c 

0.054891 




C-nat 

06000.60c 

0.737823 

AS4/3502 



H-l 

01001.60c 

0.214368 

graphite/epoxy 

1.578 

1.901141 

0-16 

N-14 

08016.60c 

07014.60c 

0.022268 

0.020364 

composite 



S-nat 

16000.60c 

0.005102 




N-15 

07015.60c 

0.000076 

Polyethylene 

0.91 

3.296703 

H-l 

C-nat 

01001.60c 

06000.60c 

0.667954 

0.332046 

B 10 -loaded 



H-l 

01001.60c 

0.624741 

0.92 

3.26087 

C-nat 

06000.60c 

0.310296 

polyethylene 



B-10 

05010.60c 

0.064963 

Be' 1 - loaded 



H-l 

01001.60c 

0.624741 

0.92 

3.26087 

C-nat 

06000.60c 

0.310296 

polyethylene 



Be-9 

04009.50c 

0.064963 




H-l 

01001.60c 

0.586502 

LiF-loaded 



C-nat 

06000.60c 

0.290897 

1.13 

2.654867 

F-19 

09019.60c 

0.061286 

polyethylene 



Li-7 

03007.60c 

0.056203 




Li-6 

03006.60c 

0.005112 




H-l 

01001.60c 

0.533709 

Li 2 CO (-loaded 



C-nat 

06000.60c 

0.298219 


1.06 

2.830189 

0-16 

08016.60c 

0.101086 

polyethylene 



Li-7 

03007.60c 

0.061392 




Li-6 

03006.60c 

0.005584 

6 



H-l 

01001.60c 

0.527656 

Li 2 C0 3 -loaded 

1.08 

2.777778 

C-nat 

06000.60c 

0.295406 

polyethylene 

0-16 

08016.60c 

0.100132 




Li-6 

03006.60c 

0.076806 


Table 2. Energy Integrated Transmitted and Reflected Neutron Flux From Each Material per Source Neutron 


Material ID 

Transmitted flux 

Reflected flux 

neutrons/cm 2 -sec -source 

neutrons/cm 2 -sec -source 

Source 

4.302 x 10" 2 


2024 A1 alloy 

1.904 x 10" 2 

4.946 x 10“ 2 

Ti6A14V alloy 

2.249 x HL 2 

4.848 x 10“ 2 

AS4/3502 composite 

1.381 x 10" 2 

5.182 x 10“ 2 

Polyethylene 

0.962 x 10" 2 

5.361 x 10“ 2 

B 10 PE 

0.850 x 10" 2 

5.187 x 10“ 2 

Be 9 PE 

0.994 x 10" 2 

5.357 x 10“ 2 

LiF PE 

1.022 x 10" 2 

5.285 x 10“ 2 

Li 2 C0 3 PE 

1.088 x 10" 2 

5.270 x 10“ 2 

Li 2 C0 3 PE 

1.041 x 10" 2 

5.221 x 10“ 2 
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Figure 3. Normalized neutron energy spectrum for transmission flux for four reference materials and for various lithium- 
containing materials. 
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